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Abstract 

The paper is devoted to the local classification of generic control-affine systems on an 
n-dimensional manifold with scalar input for any n > 4 or with two inputs for n = 4 
and n = 5, up to state-feedback transformations, preserving the affine structure (in C°° 
category for n = 4 and category for n > 5). First using the Poincare series of moduli 
numbers we introduce the intrinsic numbers of functional moduli of each prescribed number 
of variables on which a classification problem depends. In order to classify affine systems 
with scalar input we associate with such a system the canonical frame by normalizing some 
structural functions in a commutative relation of the vector fields, which define our control 
system. Then, using this canonical frame, we introduce the canonical coordinates and find 
a complete system of state-feedback invariants of the system. It also gives automatically 
the micro-local (i.e. local in state-input space) classification of the generic non-affme n- 
dimensional control system with scalar input for n > 3 (in C°° category for n = 3 and in 
C w category for n > 4) . Further we show how the problem of feedback-equivalence of affine 
systems with two-dimensional input in state space of dimensions 4 and 5 can be reduced to 
the same problem for affine systems with scalar input. In order to make this reduction we 
distinguish the subsystem of our control system, consisting of the directions of all extremals 
in dimension 4 and all abnormal extremals in dimension 5 of the time optimal problem, 
defined by the original control system. In each classification problem under consideration 
we find the intrinsic numbers of functional moduli of each prescribed number of variables 
according to its Poincare series. 
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1 Introduction 



In the paper for the convenience of the presentation all objects are C°° without special 
mentioning, although all constructions and some statements remain valid in an obvious way also 
in C k category for an appropriate finite k. On the other hand, some of the statements are known 
to be true only in the real analytic category, which will be indicated explicitly. 
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Let M be an n-dimensional manifold, fo> fi, ■ ■ ■ , fr be vector fields on M, r < n. Consider 
the following control-affine system with r inputs on M 



r 



<1 = foil) + ^u r f r (q), q G M, u\ . . . u r £ R. 



(1.1) 



k=l 



We say that a system of the type IJ1-1J) is an (r, n) control-affine system. We also assume that 
at given point qo 



Consider the group FB qo of state- feedback transformations, preserving an affine structure and 
the point q$, i.e. transformations of the type 



where $ is a diffeomorphism in a neighborhood of qo, A(q) G W, B{q) is a r x r-matrix, 
det£>(go) 7^ 0. This group of transformations acts naturally on the set of germs of systems of 
the type (|1.1|) and defines the equivalence relation, called state-feedback equivalence. The natural 
question is when two germs of the systems of the type (|l.l(l are state- feedback equivalent. 

Let us roughly estimate the "number of parameters" in the considered equivalence problem. 
The set of r-dimensional affine subspaces in R n forms (r + l)(n — r)-dimensional manifold. 
Therefore, if the coordinates on M are fixed, then the control system of the type (|l,lj) can 
be defined by {r + \){n — r) functions of n variables. The group of the coordinate changes 
on M is parameterized by n functions of n variables. So, by a coordinate change one can 
"normalize", in general, only n functions among those (r + l)(n — r) functions, defining our 
control system. Thus we may expect that the set of orbits of generic germs of systems (jl.l|) at 
qo G M w.r.t. the action of the group of transformations of the type (jl.3j) can be parameterized 
by (r + l)(n — r) —n = r{n — r — 1) arbitrary germs of functions of n variables and a number of 
germs of functions, depending on less than n variables (see the next section for the discussion 
about the number of these additional functional invariants). 

According to the last estimate the only case when the functional parameters (also called 
functional moduli) are not expected in the parameterization of the set of orbits of generic 
germs of systems IJ1-1JI is the case r = n — 1 (i.e., corank 1 control-affine systems). In this 
case under assumption (|1.2|) there is the natural one-to-one correspondence between the set of 
control-affine systems, up to feedback transformations, and the set of differential 1-forms on 
the ambient manifold: to any affine system (|1.1|) one can assign a unique differential 1-form u 
such that u(fi) = for i = l,...n — 1 and w(/o) = 1- So, the state-feedback classification 
of corank 1 control-affine systems satisfying (|1.2j) is equivalent to the well-known classification 
of differential 1-forms w.r.t. the action of the group of diffeomorphisms (see, for example, [5], 
section 3 and Appendix C there ). In particular, all germs of (|l.lj) such that the underlying 
vector distribution span(/i, . . . , f n -i) is contact for odd n or quasi-contact for even n (which is 
generic assumption) are state-feedback equivalent to the control-affine system, corresponding to 
the classical Darboux model (note also that in [Hj normal forms for codimension 1 singularities 
are given too). In the context of corank 1 control-affine systems it is not worse also to mention 
the works [3], [S]) [0]] where the case f(qo) G span(/i(go), • • • , /n-i(<Zo)) was treated. 

The case r = 1, n = 3 (the smallest dimensions, when the functional parameters appear) 
was treated in (section 3, Proposition 3.2 there). In particular, it was shown that the set of 




(1.2) 




(1.3) 
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orbits of generic germs of the systems w.r.t. the action of the group of transformations of the 
type (|1.3|) can be parametrized by one arbitrary function of 3 variables, 2 arbitrary functions of 
two variable , and the discrete invariant from the set {—1,1}. 

Remark 1 Actually, in Proposition 3.2 of the two functions of two variables satisfy 
certain conditions on coordinate subspaces of some special coordinates, which are canonical up to 
some reflections, but instead of these functions one can take their appropriate partial derivatives, 
which are already arbitrary. The functions of the parameterization are state- feedback invariants 
up to some reflections in the coordinates. □ 

In the present paper we make a classification of generic germs of systems of the type (|1.1|) . 
up to state-feedback equivalence, in the following cases 

1. r = 1, n = 4; 

2. r = 1, n > 5 in the real analytic category; 

3. r = 2, n = 4; 

4. r = 2, n = 5 in the real analytic category. 

In general, statements of the kind "the classification problem depends on the tuple of func- 
tional invariants, consisting of certain number of functions of each number of variables" need to 
be clarified: although the number of functional invariants, depending on the maximal number of 
variables, could be found in the more or less rigorous way by counting of dimensions arguments 
(as was done above), it seems that the numbers of functions of each number of variables less 
than the maximal one depend on the way of normalization of the objects in the considered 
classification problem. Besides, these numbers could be changed rather arbitrary by mixing, 
combining or separating the formal Taylor series of this functional invariants without losing any 
information (at least if we work in the category of formal Taylor series or in the real analytic 
category). In [2j (section 1 there) it was proposed to use the so-called Poincare series of the 
moduli numbers of the classification problem in order to determine intrinsically the number of 
functions of each prescribed number of variables, on which some classification problem depends. 
In section 2 below, using the Poincare series, we give a canonical selection of these numbers. The 
way the canonical parameterization is obtained indicates the presence of an interesting algebraic 
structure on the set of all tuples of fundamental invariants parameterizing given classification 
problem. For the moment, this algebraic structure remains hidden and needs further research. 

In the case of scalar input our method of the classification is similar to the procedure, used 
in pQ for the case n = 3 and it is described in section 3. It consists basically of the following 
two steps: first for any control system, satisfying some genericity assumptions, we construct the 
canonical frame by normalizing some structural functions in a commutative relation of the vector 
fields, which define our control-affine systems; then, using this canonical frame, we introduce the 
canonical coordinates and find the complete system of state-feedback invariants of the system. 

Besides, to any control system 



on an m-dimensional manifold S (the state space) with one-dimensional control space V one 
can assign the following control-affine system on the (m + l)-dimensional state-space S x V. 





(1.4) 




(1.5) 
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(here we look on u as on a new state variable, v is the new control, fo = (J-(y,v),0) T and 
f\ = (0, 1) T in the notations of It turns out that having the local classification of 

generic m + 1-dimensional control-affine systems with scalar input, one gets also the micro- local 
(i.e. local in state-input space) classification of the generic m-dimensional control system (see, 
Remark |H1 at the end of section 3). 

Further, in section 4, we show that the problem of state-feedback classification of the control- 
affine systems with two-dimensional input in dimensions 4 and 5 can be reduced to the previous 
problem for the control-affine systems with scalar input in the same dimensions. In order to 
make this reduction we distinguish the subsystem, consisting of the directions of all extremals 
in dimension 4 and all abnormal extremals in dimension 5 of the time optimal problem, defined 
by the original control system. 

In each classification problem under consideration we find the intrinsic numbers of functional 
moduli of each prescribed number of variables according to its Poincare series. 

Finally note that the problem, considered here, is different from one, considered in the paper 

which has the similar title. In the mentioned paper the authors study germs of n-dimensional 
control-affine system with scalar input at an equilibrium point qo, i.e. when fo(qo) € {R/i^o)}- 
Their method is a generalization of technique, developed in [5] and [1], which is similar to 
classical Poincare-Dulac procedure for normalization of vector fields near a stationary point. 
Therefore in the method of [7j it is crucial that qo is an equilibrium point. Here we classify the 
control-affine systems with scalar control near non-equilibrium point, which seems very natural 
in view of the fact that this classification, except the case n = 2, a priori contains functional 
moduli. Besides, the feedback invariants, constructed here for generic germs could be used also 
for the problem with equilibrium points by passing to the limit. 

2 Poincare series and the intrinsic number of functional invari- 
ants. 

We start with some terminology. Let M be a smooth manifold. Fix a point qo G M . 
Consider a set O of germs at qo of smooth objects on M (for example, Riemannian metrics, 
vector distributions, control-affine systems) such that the group of local diffeomorphisms Diff 9o , 
preserving the point qo, acts naturally on it. This action defines the equivalence relation on O. 

Denote by J k {0) the space of all fc-jets at qo of objects from the set O. We say that the 
set O C O is a generic subset of O if there exists an integer k > and a Zariski open set U in 
J k (0) such that 

= {beO : jfc-jet of b belongs to U}. 

By classification problem on O we mean the problem to find a system of fundamental in- 
variants for objects from some generic subset of O such that two generic objects are equivalent 
if and only if they have the same systems of fundamental invariants. Let O be a generic subset 
of O, which is invariant w.r.t. the action of the group Diffgy. 

Definition 1 A mapping F from the set O to the set C£°(1R Z ,M) of germs at of smooth 
functions in W, which is invariant w.r.t. the action of the group Diff go on O, is called a 
functional invariant of I variables of a generic subset of objects from O. 

When the object b £ O is fixed, we will mean by the functional invariant also the value of the 
mapping F at b. We will denote this germ of function by the same letter F without special 
mentioning. 



4 



Now let us describe the Poincare series of the moduli numbers of the classification problem. 
The action of the group Diff 5o induces the action Ak of some finite dimensional Lie group Gj~ 
on the space J k {0) for any integer k > 0. So , Ak is a mapping from Gk x J k (0) to J k (0). 
Given any b € J fc (C) let „4f be the mapping from G k to J k (0) such that ,4|(-) = A k (; b). Let 
efc be the identity of the group Gk- Set 

m(k) = dim J k ((D) — max rankdAufejfc) = min corank dAu(ek) (2-1) 
be.J k (o) be.J k (0) 

(here rank(L4|!(efc) and corank dA\{ek) are the rank and the corank of the differential of A\ 
at e& respectively). Roughly speaking m{k) is the dimension of the space of orbits w.r.t. the 
last action Ak- The number m(k) is called the moduli number of the k-jets. The Poincare 
series of the moduli numbers of the classification problem (or shortly the Poincare series of the 
classification problem) is by definition the following function: 

oo 

M(t) = ^2m(k)t k (2.2) 

fc=0 



Remark 2 Since the integer- valued function b i— > ranke&Auejfe) takes its maximal value at 
a Zarisky open set, in (|2.1|) we can replace O by any its generic subset O. □ 

The Poincare series could be useful in evaluating of the number of the functional invariants 
of the given number of variables, on which the given classification problem depends, because of 
the following well-known fact: if one denotes by ji(k) the dimension of the space J fc (M',IR) of 
fc-jets of functions of /-variables, then the corresponding Poincare series of numbers ji{k) satisfies 

k=o v ; 

(here one uses that ji{k) = ;Tfc! )■ So, if , for example, the Poincare series of some classification 
problem is equal to 

n 

1=1 v ; 

where all pi are nonnegative integers, then it is natural to conclude that this problem depends 
on the tuple, consisting of pi functional invariants of % variables for each 1 < % < n, while the 
parameter w (i.e. the order of zero of the Poincare series M(t) at t = 0) is equal to the minimal 
k > such that the action of the group Gk on the space of J k (0) is non-transitive. 

Till now nothing is known about the form of the functions M(t) for general classification 
problems. For example, the following open question is stated in [2j: Is it true that the Poincare 
series of moduli numbers are rational functions in the most classification problems? In the next 
sections we will show by direct computations that in all classification problems 1-4 listed in 
Introduction it is true and moreover the function M(t) has a unique pole at t = 1. On the other 
hand, in all considered cases (except the cases r = 1, n = 3 or 5) the Poincare series has no a 
representation of the type Q2.4JI with nonnegative pi- Below we give an algorithm to extract the 
number of functional invariants from Poincare series also in these cases. 

From now on we will suppose that the Poincare series M(t) of the classification problem is a 
rational function with a unique pole at t = 1. Let wq be the order of zero of the function M(t) 
at t = 0. 
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Lemma 1 For any integers w >wq and I there exist a unique polynomial R(t) with 

deg R(t) <w-w (2.5) 
and a unique rational function Q(t) with the unique pole at t = 1 such that 

t w °R(t) 

M(t)= {1 _ t yi 1 +t W Q(t)- (2-6) 

Proof. Let us fix I £ Z and prove the existence of a representation of the type (|2.6|) for any 
w > w o by induction in w. 

If w = wo, then from the condition (|2,5|) it follows that R(t) = 0. Then by definition of order 
of zero the function Q(t) = -^M{t) is rational with the unique pole at t = 1, which implies 
(ESJ). 

Now suppose that a representation of the type (|2.6j) exists for some w = w, w > wq, and 
prove its existence for w = w + 1. For this let Q(t) and R(t) be as in the representation (|2.6I) 
for w = w. Denote by 

t r ltJ (i-t)' 



Qi(t) = 7 (Q(t)- ,, V t+1 . (2.7) 



Then by construction Qi is also the rational function with the unique pole at t = 1. Expressing 
Q(t) from (|2.7|) and substituting it into (|2.6|) . one has 

t W0 (R(t) +Q(0)t™- wo ) 
M (*) = (i_t)t+i - + (2-8) 

Since deg(i?(t) + Q(0)t M '~ M ' ) < — wo + 1 it implies the existence of a representation (j2.6|) also 
for it; = tu + 1. This completes the proof by induction of the existence part of the lemma. 

Now let us prove the uniqueness part. If there exists another representation of M{t) of the 
type ()2.6(1 with a polynomial R(t), deg R(t) < w — wo, and a rational function Q(t) instead of 
R(t) and Q(t), then we have the following identity 

l+i 



R(t) - R(t) = t w - wo (l - t) l+1 (Q(t) - Q(t)) . 

It implies that the polynomial R(t) — R(t) has zero of order not less than w — wq. On the other 
hand, by assumptions deg(i?(t) — -R(i)) < w — wq, which implies that R(t) = R(t) and then also 
Q{t) = Q(t). □ 

We will call the representation (|2.6j) (with R(t) satisfying 1)2.51) ) the (w,l)— representation 
of the function M(t). Let N be the order of pole of (1 - t)M(t) at t = 1. 



Definition 2 The (w , I) -representation \2.b}) of M(t) with R(t) and Q(t) satisfying 

w—wq — I N 

£ nt\ Q(t) = J2 

i=0 j=h 

is called nice, if 1 < I < N , l\ > I, and all coefficients ri, qj in (|2.9|) are nonnegative integers. 



R(t)= ^ n t\ Q(t) = g 9Jl ^0 (2. 
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Of course, in general a rational function M(t) = rfz^w+i > where is a polynomial (even 
with integer coefficients), may not have any nice (w, ^—representation. But if the function 
M{t) is the Poincare series of a classification problem, which can be parameterized by functional 
invariants in some reasonable way, then M{t) has at least one nice representation. 

To be more precise and to explain why the nice representation of the Poincare series are 
interesting let us introduce an additional terminology. Let F be a functional invariant of I 
variables of a generic subset O of objects from O. Denote by irk the natural projection from the 
set O to the space J k (0). Let (x\, . . . xi) be the standard coordinates in M. 1 . Let Co°(R z ,M) be 
the set of germs at of smooth functions in M 1 , Given a tuple . . . , i{) of nonnegative integers 
denote by the mapping from Cjf(R l ,R) to R, which assigns to any function / :! h| 
the value of its partial derivative — {- at 0. 

Definition 3 We say that the functional invariant F of I variables of a generic subset O C O 

has the weight w if for any integer k > w the following condition holds: for any tuple (ii, . . . , i{) 

i 

of nonnegative integers with y^i s < k — w there exists a mapping ©(^...j;) : J k (0) i— > R such 

8=1 

that the following diagram 

O £-C§°(K',K) 



I 

is commutative, while for any tuple of nonnegative integers with > k — w such 

s=l 

mapping <&i lr ..i t does not exist. 

Essentially the weight of the functional invariant F is the integer w > such that for any k > w 
the fc-jets of our objects depend on all derivatives of F of order not greater than k — w but do 
not depend on derivatives of F of order greater than k — w. 

The classification problem with functional moduli will be called regular, if there exist integers 
w\ and W2, n%, and n2 ( where < w\ < W2, 1 < n\ < n2) and a {ri2 — n\ + 1) x (u>2 — w\ + 1) 
matrix P with nonnegative entries pij such that the set of orbits of generic germs of systems 
w.r.t. the group Diff^, can be parameterized by the set of functional invariants , consisting 
of pij functional invariants of i + n\ — 1 variables and weight j + ui\ — l for any 1 < i < n2 — n\ + \ 
and 1 < j < W2 — W\ + 1. Without loss of generality we suppose also that on the first and the last 
column and row of the matrix P there is at least one nonzero entry. Directly from the definition 
of weight, formula (|2.Hj) . and Remark El it follows that the Poincare series of such classification 
problem satisfies 

102—101+1 /"a— ii+l \ 

M(t) = / - y ip y - %- . (2.io) 

Note that w\ = wq and ri2 = N, where as before wq is the order of zero of M{t) at t = and A'" 
is the order of pole of (1 — t)M(t) at t = 1, but all other parameters, appearing in 1)2. lOj) could 
not be uniquely recovered from the Poincare series M(t). In the considered situation we will say 
that the classification problem admits a regular {w2,n\) -parameterization with parameterization 
matrix P. 



J k (0) 
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Given some (1V2, ni)-parameterization of the classification problem one can easily build new 
(1V2, ni)-parameterization with another parameterization matrices. Indeed, take some functional 
invariant F of the weight jo, depending on iq variables, say x%, X2, ■ ■ ■ , Xi , where 2 < iq < 
112 — n\ + 1 and 1 < jo < 102 — w\. Let G{x\, . . . , X{ ) be the function such that 

F(xi, . . .,x io ) = F(xi, . . . ,x io -i,0) + x io G(x 1} . . .,x io ). (2.11) 

Then we can obtain the new parameterization of the classification problem by replacing the func- 
tional invariant F(x\, . . . ,Xi ) by two functional invariants F(x\, . . . , Xj _i, 0) and G(x\, . . . , Xi ). 
Obviously, the first invariant has weight jo and depends on iq — 1 variables, while the second 
one has weight jo + 1 and depends on iq invariants. The matrix of the new parameterization is 
obtained from the original one by decreasing the (iq, jo)-entry by 1 and increasing both (io~ 1, jo)- 
entry and (io, jo + 1) entry by 1. Such transformation on the set of (N — m + 1) x (1U2 — u>o + 1) 
matrices will be called an elementary transformation. Conversely, given two functional invari- 
ants G\ and G2 such that G\ depends on iq — 1 variables say x±,... ,Xj _i and has the weight 
jo, while G2 depends on iq variables say xi,... ,Xi and has the weight jo + 1 (here again 
2 < «o < n 2 — n\ + 1 and 1 < jo < W2 — wi) one can build the new parameterization by replacing 
the invariants Gi and G2 by one invariant G\ + a^G^, which depends on iq variables and has 
the weight jo. Of course, in this case the matrix of the new parameterization is obtain from the 
original one by the transformation, which is inverse to the elementary one. 
Now for convenience denote 

Ki = N-m + l, K 2 = w 2 -w + l. (2.12) 

Note that among all matrices, which can be obtained from the given K\ x K2 matrix P with 
integer entries by a composition of a finite number of elementary transformations and their 
inverses, there exists a unique matrix, denoted by Norm(P), such that all its entries, except those 
lying on the first row and the last column, are equal to zero. To prove the existence of Norm(P) 
one can vanish the entries of the matrix P by a composition of elementary transformations and 
their inverses step by step, starting from the entry in the left-lower corner, going along the first 
column from the bottom to the top till the entry on the second row , then passing to the bottom 
of the second column, going along it from the bottom to the top till the entry on the second row 
and so on till the column before the last one. The uniqueness follows from the fact that if we 
put the entries of the matrix Norm(P) instead of the entries of P in the representation (|2.10[) . 
then we obtain the (w2, ni ^representation of the Poincare function M(t). This representation 
is unique according to Lemma ^ and the matrix Norm(P) is obviously uniquely recovered from 
it. Also it is not difficult to express all nontrivial entries of Norm(P) by the entries of P: 

(Norm(P)) w = + £ £ ( , jPk+lj-l, l<j<K 2 -l, (2.13a) 

1=0 k=i ^ ' 

Ki—i+l K 2 -l sj^ _ ^ j _ -|\ 

(Norm(P))^ 2 = Pi , K2 + E ( 2 / W<>fc> 2<i<Ki, (2.13b) 

1=0 k=i ^ ' 
(Novm(P)\ K2 =p ljK2 (2.13c) 

The last relation can be proved, for example, using the procedure of passing from P to Norm(P), 
described above, and the following well-known combinatorial identity: 

/ i + k — 1\ fn + k 



E 

i=i 



fe + 1 
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If the matrix P has only nonnegative integer entries, then the matrix Norm(P) is obtained 
from P by a finite composition of elementary transformations (without using their inverses) and 
has also only nonnegative integer entries (which follows also from relations l|2.1H|) ). Moreover, if 
we put the entries of the matrix Norm(P) instead of the entries of P in the representation ()2.1(jj) . 
then we obtain the nice (w2, ni ^representation of the Poincare function M(t) of our classification 
problem. We also say that this nice representation corresponds to the matrix Norm(P). We can 
summarize all above in the following 

Proposition 1 If the classification problem admits a regular (w2,ni) -parameterization with 
parameterization matrix P, then it admits a regular (w^n\) -parameterization with parame- 
terization matrix Norm(P) and its Poincare series has the nice (w2,n\) -representation, which 
corresponds to the matrix Norm(P) . 

The last proposition indicates that the nice representations of the Poincare series (if they 
exist) may be used in the definition of the intrinsic number of functional invariants of each 
number of variables and weight, on which the given classification problem depends. Suppose 
that the Poincare series M(t) has the nice representation for some (w,l). So, the set 

NS(M(t)) d = {(w, I) : the (w, I) -representation of M(t) is nice} (2.14) 

is not empty. The natural question is what pair to choose from NS (M(t))? To answer this 
question we propose to introduce the order -< on the set of ordered pairs (w, I) in the following 
way: (w, I) -< (w, I) if and only if w < w or w = w, but I > I. By Definitional 

NS(M(t)) C {(w, l):w>w ,l< N}, (2.15) 

which implies immediately that the set NS(M(t)) contains the minimal element w.r.t. the intro- 
duced order -<. This minimal element will be called the characteristic pair of the classification 
problem. Denote it by (w,l). Let C be the (N — I + 1) x (w — wo + 1) matrix such that the 
(w, ^-representation of M{t) corresponds to the matrix C. 

Definition 4 The (ij) -entry of the matrix C is called the intrinsic number of the functional 
invariants of i + 1 — 1 variables and the weight j + WQ — 1 of the considered classification problem. 
The matrix C is called the characteristic matrix of the classification problem. Any regular (w, l)- 
parameterization of the problem (if it exists) with the parameterization matrix C is called the 
characteristic regular parameterization. 

In general, it is better to have a parameterization, consisting of invariants, which have mini- 
mal possible weight and depend on maximal possible number of variables. Our definition of 
characteristic parameterization is in accordance with this goal. Actually the maximal weight of 
invariants, appearing in characteristic regular parameterization is not greater than the maximal 
weight of invariants, appearing in any other regular parameterization. Besides the minimal num- 
ber of variables in invariants of regular parameterization is not less than the minimal number of 
variables in invariants of any other regular parameterization, having the same maximal weight 
of invariants characteristic one. 

One can improve the formula (|2.15j) for the localization of the set NS(M(t)). Indeed let d 

be the degree of the rational function M{t) at infinity. Namely, if M{t) = > where Qi(t) 
and Q2(t) are polynomials, then d = degQi(i) — degQ2(t). Then 

NS(M (t)) C {(w, I) : w > wo, 1 < I < min (w - d - 1, N)}. (2.16) 
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To prove (|2.16|) we actually have to prove that if the pair (w, I) € NS(M(t)) then I < w — d — 1. 
Indeed, from Q2.6|) and Q2.9|) it follows that d = max (w — — 1, i^o + deg R — l\ — 1). But first, 
since Zi > I, we have to — l\ — 1 < it; — / — 1 and secondly, since degi? < w — wo, we have 
u> + deg R — l\ — 1 < w - I - 1. Therefore d < w — I — 1, Q.E.D. 
From (|2.16j) it follows also that 

NS(M(i)) C {(iM) : w > max(w ,d + 2)}. (2.17) 

The relations (|2.16|) and (|2.17j) may be useful in searching for the characteristic pair of the 
classification problem. Another useful property of the set NS(M(t)) can be formulated as 
follows: 

Lemma 2 Assume that the function M (t) has the nice (w, I) -representation \2. fij) . the func- 
tions R(t), Q(t), and the number l\ are as in 12. y\) . and l\ = I (or, equivalent, qi > 0), then 
(w- 1,1-1) 0NS(M(t)). 

Proof. Let S(t) be the polynomial such that M(t) = r^zqw+i ■ Then, using the assumption 
I = 1%, it is easy to get 

deg S{t) =w + N-l. (2.18) 
Moreover, directly from ()2.6|) and (|2.9I) one can obtain that 

1¥ ^ r[ =(-l) N -\ W + N-l)!q l . (2.19) 
On the other hand, if the (w — 1,1 — l)-representation of M(t) has the form 

then deg S(t) = max(tu — 1 + iV — l2,w + N — I — 1). Comparing this with (|2.18j) one gets easily 
that I2 = I — 1- But then by analogy with ()2.19j) (applied for (w — 1, 1 — l)-representation instead 
of (w, /)-representation) one has 

d^rr = (-if- l+ \™ + N-t)\ (2.21) 

Comparing (|2.19|) and (|2.21|) . we obtain that qi—\ = —qi- Hence qi-\ < and the (w — 1,1 — 1)- 
representation (|2.20|) is not nice. □ 

As a direct consequence of Proposition ^ the previous lemma, and the relation (|2.13c|) one 
has the following 

Corollary 1 // the classification problem with the Poincare series M{t) admits the regular 
(w, I) -parameterization with the parameterization matrix P such that the entry in the right- 
upper corner of P is positive (i.e., in the previous notations pi )W - WQ +i > 0), then (w — 1,1 — 1) 
NS(M(t)). 

In the sequel we will show that all classification problems 1-4 listed in Introduction are regu- 
lar. For each of this problems we will describe explicitly some its regular (u>, Z)-parameterization 
such that (w, I) is the characteristic pair of the problem and find its characteristic matrix, which 
also gives the way to obtain the characteristic parameterization. 
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Remark 3 Sometimes (as in the case of (1,3) control-affine systems) in the parameterization 
of a classification problem appear invariants from the finite set in addition to the functional 
parameters. Also these functional parameters may be not invariants by themselves, but up 
to some finite group of transformation of coordinates. We call such classification problems 
quasiregular. Obviously, the additional invariants from the finite set and the additional freedom 
in choosing the coordinates, up to the finite group of transformation, do not affect on the Poincare 
series. Therefore the identity (|2.1fl|) for the Poincare series is true also for quasiregular problems, 
where, as before, pij is the number of functional invariants of i + m — 1 variables and weight 
j+wi — 1 for any 1 < i < 712— ni+1 and 1 < j < vfy—wi+l, and we can extend to the quasiregular 
case the same terminology and the constructions, as in the regular one. In particular, according 
to Proposition 3.2 of (see also Remark ^ and the paragraph before it in Introduction), the 
state-feedback classification problem for (1,3) control-affine system admits quasiregular (2,2)- 

parameterization with the following 2x2 parameterization matrix P = ^ ^\ . Besides, wq = 1 

and N = 3. By the inverse to the elementary transformation one can transform P to the matrix 



It is not difficult to see that by erasing the last column of the matrix P one obtains the char- 
acteristic matrix C = (1, 1) T of the considered classification problem and the characteristic pair 
is equal to (1,2). So, the characteristic quasiregular parameterization consists of one function 
of 3 variables and the weight 1, one function of 2 variables and the weight 1, and the discrete 
invariant from the set {—1,1}. This parameterization is obtained from the original one by a 
rearrangement of the invariants, which corresponds to the inverse to elementary transforma- 
tion, transforming the matrix P to the matrix P (such rearrangements were described in the 
paragraph after the formula (|2.11[) ).D 

3 Classification of (l,n) control-affine systems for n > 4 

For r = 1 the system p. 1(1 has the following form 





? = /o(g)+«/i(g), ?6M, aei. 



(3.1) 



Our genericity assumptions are 



dimspan(/ , fx, /o], • • • (ad/i) n 2 / ) = n, 
dimspan(/i, [/ , [/ , fx]] , [fx, fo], ■ ■ ■ (ad/i) n ~ 2 /o) = n, 



(3.2) 
(3.3) 



and for n > 5 also 



dimspan(/ , /i, [/ , [fo, fx}] , [fxJo], ■ ■ ■ (ad/i) n 3 / ) = n. 



(3.4) 



The group of feedback transformations 



u = P(q)u + a(q),p(q) ^0 



(3.5) 
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acts naturally on the set of pairs of vector fields (/o,/i). The orbit w.r.t. this action is 

O {fojl) = {(/o + aft, @h) : a, (3 : M ^M. are functions, (3 / 0}. (3.6) 

The first observation is given by the following 

Proposition 2 If the pair (/o,/i) satisfies conditions KH. ty) and iV..V)) . then there exists a 
unique pair (Fq,Fi) G Q(f f-A such that 

n-2 

[Fq, [F , F 1 )) =F + hF L + I 2 [F 1} F ] + ^ I^dF^Fo (3.7) 

k=3 

Proof. By assumptions ()3.2j) the vector fields fo, fx, [fi, fo], ■ ■ ■ (ad/i) n_2 /o constitute the 
frame on M. Therefore there are functions E, N, J\,..., J n -2 such that 

n-2 

[/o,[/o,/i]] =Ef + Jxh + J2[fi,M + Z[f 1 ,[fi,M] + ^J fc (ad/i) fe / ( 3 - 8 ) 

fc=3 

Take some pair (/ ,/i) G 0^ f()Jl} , 

/o = /o + «/i, /i = /3/i. (3.9) 

Suppose that 

n-2 

[fo, [fo, fi}} = Efo + Ji/i + MfiJo] + Z[fi, [A, fo]] + U^hfh (3.10) 

fc=3 

First note that 

E = (3E. (3.11) 
It follows immediately from ()3.9j) and the following relations 

[fo, [/o,/l]]] = P[fo, [/o,/l]] - aPfa, [fi,fo]] modspan^x, [f ,fi]), (3.12) 

(ad/i)*/o G span(/i, . . . , (ad/a) fc / ), fc G N. (3.13) 
From assumption ()3.3j) it follows that E ^ 0. Therefore, taking (3 = -g, we make 

£ = 1. (3.14) 

Let us denote by O(f 0> f{) the set of all pairs (fo, fi), satisfying (|3.14l) . We can assume from the 
beginning that the original pair (fo, f\) belongs to O^ j 1 - ) , i.e. E = 1 (we make this assumption 
just in order to avoid extra notations). If (fo,fi) G Off Q j x -\, then also E = 1. Hence (3 = 1 
or equivalently f\ = f\. In other words, condition ()3.14|) normalizes the vector field f\ or the 
direction defining the straight line in the set of admissible velocities of the system ()3.1JI at any 
point. 

Further, from 1)3.9(1 . taking into account that = 1, it follows easily that 
(ad/i)*/o = (ad/!) fc /o modspan(/ 1 ), k G N. 
This and relation ()3.12() imply that 

Z = Z-a. (3.15) 

Setting a = Z, we make Z = 0, which normalizes the drift fo. So, we have proved that there 
is a unique (fo,fi) G O(f j{\ such that E 1 = 1 and Z = 0, which completes the proof of the 
proposition. □ 
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Remark 4 The mappings h, ■ ■ ■ ,I n ~2 from M to R, denned by identity ((3.7|) . are state- 
feedback invariants of the control system ((3.1() . □ 

The vector field Fq and the pair of vector fields (.Fo > Fi ) from Proposition [21 are called i/ie 
canonical drift and i/ie canonical pair of the system 1)3.1(1 respectively. 

Remark 5 Actually, in the case n = 4, the vector Fo(q) is the velocity of the unique 
abnormal extremal starting at q of the time optimal problem defined by system 1)3.1(1 . □ 

Now fix some point go £ M. Denote by e l f the flow, generated by the vector field / and 
q o e 1 ^ the image of the point q w.r.t. this flow. Let <& n : R n HMbe the following mapping 

<S>4{xi,x 2 ,x 3 ,x 4 ) = q o e '*[ F ^M o e^i^o] g*a*i e ^ 0) (3 16) 



<M*1, . . . , x n ) = q o e^F^Fo o . . . e *5 K^o]] ^ 

From assumption ((3.2(1 in the case n = 4 or assumption ((3.4(1 in the case n > 5 it follows 
that 3>^(0) is bijective. Hence <&~ l defines the canonical coordinates in a neighborhood of qo (or 
shortly the canonical coordinates at qo). Denote 

2- fc = I fc o* n , fc = l,...n-2. (3.18) 

Assigning to any generic germ at go of control-affine systems ((3.1(1 the function we obtain 
the functional invariant of n variables of this set of objects in the sense of Definition ^ for any 
1 < k < n - 2. 

Now let us consider the cases n = 4 and n > 5 separately: 

a) The case n = 4. By (|3.16() and ((3 .7(1 . in the canonical coordinates the vector fields Fo 
and F\ have the following form: 

fc=i 

where the components of i*\ satisfy the following second order linear ordinary differential equa- 
tions w.r.t. the variable x\ 

^+l 2 ^±-l iak -5 hk = k = 1,2,3,4, (3.20) 
oxf OXi 

with the following restrictions on the initial conditions for any k = 1,2,3,4 

a k (0, x 2 ,x 3 , x 4 ) = 5 2k , (3.21a) 
dak 

— (0,0,x 3 ,x 4 ) = -<5 3fc , (3.21b) 
d 2 a k 

-(0,0,0,x 4 ) = -hk, (3.21c) 



dx\dx 2 



13 



where <5jj is the Kronecker symbol. Let for any k = 1, 2, 3, 4 

/3 fc (x2,x 3 ,x 4 ) d = ^y^ (0,X2,X3,X4), (3.22a) 
d e / d 3 ak 

^k[X3,X4) = - — ^ — ^— (0,0,rc 3 ,x 4 ). (3.22b) 
0x10x20x3 

So, with any germ at qo of a four-dimensional affine system (|3.1j) . satisfying genericity assump- 
tions (|3.2|) and ()3.3|) . one can associate the ordered tuple 



(l 1 ,l 2 ,Pl,P2,P3,Pi, tpi,ip2,ip3, ^4), (3.23) 

of state-feedback functional invariants, consisting of two germs I\ and I2 of functions of four 
variables at 0, four germs Pi, P2, P3, Pi of functions of three variables at 0, and four germs 
"01 > ^2 > "03 > "04 °f functions of three variables at 0. We call it the tuple of the primary invariants 
of the (1,4) control- affine system (|3.1j) at the point go- Note that by (|3.22j) the functional 
invariants Pk and V'fc have the weight 3 for any 1 < k < 4, while by ()3.2(Jj) and (|3.21|) the 
functional invariants I\ and I2 have the weight 2. 

Further, fixing /?& and and using (|3.21b|) and (|3.21c|) . one can find §f^(0, X2, X3, X4) for any 
1 < /c < 4 by the appropriate integrations (see (|3.31|) below). If in turn one fixes also X\ and Z2, 
then from the knowledge of §f^(0, £2, 23, 24), condition (|3.21aj) . and differential equation (|3.2()j) 
we can recover the functions afc(xi, X2, X3, X4) and therefore our control-affine system itself, just 
using the standard existence and uniqueness results from the theory of ordinary differential 
equations. We summarize all above in the following: 




Theorem 1 Given two arbitrary germs T\ and X2 of functions of four variables at 0, four 
arbitrary germs Pi, P2, P3, Pi of functions of three variables at 0, and four arbitrary germs 
ipi, tp2, ip3 , "04 of functions of two variables at there exists a unique, up to state-feedback trans- 
formation of the type hl.ty) . four- dimensional control-affine system with scalar input, satisfying 
genericity assumptions and \3. °J\) . such that the tuple (li,l2, Pi, P2, P3, Pi, ipi, 1P2, ip3, ipi) is 
its tuple of the primary invariants at the given point go ■ In other words, the tuples of the primary 
invariants give the regular (3, 2) -parameterization of the considered classification problem with 
the following 3x2 parameterization matrix P: 



(3.24) 



The Poincare series M(t) of the considered classification problem satisfies 

2t 2 ( 4 4 \ 

= + ' s [ jr^tf + jt^w) ■ (3 ' 25) 

It turns out that (3, 2) is the characteristic pair of the considered classification problem. 
Indeed, let as before wo be the order of zero of M(t) at t = 0, N is the order of pole of 
(1 - t)M{t) at t = 1, and d is the degree of M(t) (at infinity). Then from (l3~251) it follows that 
w = 2, N = 4, and d = 0. Hence from l|2~T5|) (or (j2H|)) it follows that 

NS(M(t)) C {{w, I) :w> 2}. (3.26) 
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Further, using (|2.16|) . one has 

NS(M(t)) n{(w,l) : w = 2} C {(2,1)}, NS(M(t)) D {(w, I) : iu = 3} C {(3, 1), (3, 2)}. (3.27) 

But by the previous theorem our classification problem admits (3, 2)-parameterization such that 
its parameterization metric has a positive entry in the right-upper corner. Therefore from 
Corollary [U it follows that (2,1) £ NS(M(i)). Since (3,2) -< (3,1) we can conclude from 
(f3~2l)|) and (|3~27|) that (3,2) is the minimal element of NS(M(t)). In other words, (3,2) is the 
characteristic pair of our classification problem. Also, it implies that the characteristic matrix 
C of the problem is equal to Norm(P), which can be found easily by the series of elementary 
transformations. Namely, 



C = Norm(P) = 




(3.28) 



Conclusion 1 The characteristic parameterization of (1,4) control- affine systems, up to 
state-feedback transformations, consists of 2 functional invariants of 4 variables and the weight 
3, 6 functional invariants of 3 variables and the weight 3, 2 functional invariants of 2 variables 
and the weight 2, and 4 functional invariants of 2 variables and the weight 3. 

In order to obtain a characteristic parameterization from the parameterization by the tuple of 
the primary invariants one can implement some series of rearrangement of the primary invariants 
according to the series of elementary transformations from the matrix P to Norm(P), as was 
described in the previous section (see, for example, formula Q2.11JI and the paragraph after it). 

We finish the treatment of the case of (1, 4)-affine control systems by writing the local normal 
form of such system, up to state-feedback transformation, in terms of the tuple of their primary 
invariants: Let iV be the solution of the following non-homogeneous second order linear ordinary 
differential equation w.r.t. the variable x% with prescribed initial values 



( d 2 N dN 

dx\ 2 dx\ 



X\N — 1 = 0; 
dN 



N(0,x 2 ,x 3 ,x 4: ) = 0, - — (xi,x 2 ,x 3 ,x 4 ) 

ox l 



(3.29) 



0. 



xi=Q 



and the functions p%, p 2 be the solution of the following homogeneous second order linear ordinary 
differential equations w.r.t. the variable x\ with prescribed initial values 



d 2 p t dpi 

2 +l2 ~ ^Pi = °> 



1,2, 



dx'f ' * dxi 

pi(0,X 2 ,X 3 ,X4) pi(0,X 2 ,Xs,X4) 
s | r Pl(0,X2,X3,X 4 ) ^P2(0,X 2 ,X 3 ,X4] 



1 
1 



(3.30) 



Let also 



(rx3 \ rx2 

-b~Ak + J ipk(y,X4:)dyj + J {x 2 - y)(3 k (y,x 3 ,X4,)dy (3.31) 

for 1 < k < 4 (actually Bk(x2, x 3 , x 4 ) = §^-(0, x 2 , »3, X4), where the functions are as in 
(|3.19j) ). Then a four-dimensional control- affine system ()3.1|) with the tuple of the primary 
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invariants (Ii,l 2 , Pi, P 2 , fiz, Pa, ipi, tp2, tp3, tpi) at the point qo is state-feedback equivalent to the 
following system: 

ii = l + (N + B lP2 )u 
< ± 2 = (pi + B 2 p 2 )u , «£R. (3.32) 
±i = Bip 2 u, i = 3,4 

b) The case n > 5. By (j3.17j) and (|3.7j) . in the canonical coordinates the vector fields Fq 
and F\ have the following form: 



^ ^ d ^ d 

fc=i ft 



(3.33) 



where the components a k of .Fo satisfy the following system of partial differential equations 

d 2 a k dai da k 



d n 2 a k , d^a k 



j=3 



dx[ jr{ ^ dx{ dxi dx 



+ ^2^- + a k + Iih,k = 0, 

OX i 



(3.34) 



k = 1, . . . , re, 

with the following restrictions on the boundary conditions for any 1 < k < n 

a k (0,x 2 , . . . ,x n ) = 5 2k , 
dak, 



dx\ 
d 2 a k 
dx\dx 2 



(0,0, X 3 . . . ,X n ) = 5; 



?>k ■ 



(0,0,0, £4, ...,x n ) = -5 4k , 



d 3 a k 
dx\ 



(0, . . . ,0,x j+ i, ...,x n ) = Sj +3tk , 2 < j < n - 3, 



(3.35a) 
(3.35b) 

(3.35c) 

(3.35d) 



where 6™ is the Kronecker symbol. Note also that the genericity assumption ()3.4|) implies that 



Zn-2 + 0. 

Let us introduce the following functions for any 1 < k < n: 

def d^a k 

P k (x 2 ,...,x n ) = - — — - 1 (0,x 2 , . . . ,x n ), 
OX1OX2 

def d^a k 

tp k (x 3 ,...,x n ) = - — - — — (0,0,x 3 ,... ,x n ), 



<l>kji(xi, ■■■,x n ) 



dx\dx 2 dx 3 
def d j+1 a k 



(3.36) 

(3.37a) 
(3.37b) 



dx\dxi 



a k (0,...,0, Xl ,...x n ), 2 <j < n-3, 2 < I < j + 2. (3.37c) 



So, with any germ at qo of an n-dimensional affine system (j3.1j) . satisfying genericity assumptions 
(j3.2|) and (|3.3j) . one can associate the ordered tuple 



{{ls{xi, ■ ■ -,X n )} n s= l, {Pk{x 2 , ■ ■ .,X n )} k=1 , {lpk(x 3 , ■ ■ ■ ,X n )} k=ll 

{(j) kj l{xu ■ ■ ■ , x n ) : 1 < k < n, 2 < j < n - 3, 2 < / < j + 2} 



(3.38) 
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of state-feedback invariants, of functions of n — 1 variables at 0. We call it the tuple of the 
primary invariants of the (l,n)-affine control system (|3.1[) with n > 4 at the point q^. Note 
that by ()3.22j) for any 1 < k < n the functional invariants (3k and ipk have the weight 3 , the 
functional invariants (fikji have the weight j + 1, while by (|3.20|) and ()3.21j) for any 1 < s < n — 2 
the functional invariants I s have the weight n — 2. 

Further, fixing (3k and V'fc and using ()3.35b|) and ()3.35c|) . one can find ^-(0, X2, ■ ■ ■ , x n ) 
for any 1 < k < n by the appropriate integrations. Similarly, fixing {4>kjl}{=2 ^ or gi ven 3i 
2 < j < n — 3, and using (|3.35d|) . one can find ^p^(0, X2, • • • , x n ) for any 1 < k < n by the 
appropriate integrations. Finally, if we suppose that all functions (3^, ip^, and (pkji are real 



analytic and fix also real analytic {2 S }™ = {, then from the knowledge of ^f-(0, X2, ■ ■ ■ , x n ) for all 



1 < j < n — 3, condition (|3.35a|) . and differential equation ()3.34|) we can recover the functions 
ak(x\, . . . , x n ) and therefore our affine control system itself, just using the classical Cauchy - 
Kowalewsky theorem for system (|3.34|l . We summarize all above in the following: 



Theorem 2 If n > 5, then given an arbitrary tuple ([3.38ft of real analytic functions n — 2 
germs {Ij^Zi of real analytic functions there exists a unique, up to state-feedback real analytic 
transformation of the type hl.rfy . n-dimensional real analytic control- affine system with scalar 
input, satisfying genericity assumptions ici.ty) . \S. and \3.4\ ),such that the tuple ()3.38|) is its 
tuple of the primary invariants at the given point qo . In other words, the tuples of the primary 
invariants give the regular (n — 2, 2) -parameterization of the considered classification problem in 
the real analytic category with the (n — 1) x (n — 4) parameterization matrix P such that for 
n = 5 

P= (5, 10, 10,3) T (3.39) 

and for n > 5 



On \ 




n 
2n 
2n 
V0 



n 
n 
11 
11 







n 
n 
n 
n 
ii 



(3.40) 



2/ 



(in the matrix P all entries in the triangle with vertices in the (1, l)-entry, the (l,n — 5)-entry, 
and the (n — 5, l)-entry are equal to 0, while all entries in the triangle with vertices in the 
(1, n — A)-entry, the {n — 4, l)-entry, and the {n — 4, n — 4)-entry are equal to n). The Poincare 
series M(t) of the considered classification problem satisfies 



M(t) = nt 3 



■ + 



[(l-t) n (l-t) n - 1 

n— 2 n— 1 



1 




(3.41) 



i=4 j=n—i 



(i - ty+ l (i - t) n 



+i • 



It turns out that (n — 2, 2) is the characteristic pair of the considered classification problem. 
Indeed, let as before wq be the order of zero of M(t) at t = 0, iV is the order of pole of (1— t)M(t) 
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at t = 1, and d is the degree of M(t) (at infinity). Then from (|3.25f) it follows that vjq = 3, 
N = n, and d = n — 5. Hence from (j2.17|) it follows that 

NS(M(t)) C {O,0 :w > n-3}. (3.42) 



Further, using (|2.16|) . one has 

NS(M(i)) n {(«,, : «, = 2} C {(n - 3, 1)}, 

NS(Af (t)) n {(w, l) : w = n - 2} C {(n - 2, 1), (n - 2, 2)}. 

But by the previous theorem our classification problem admits (n — 2, 2)-parameterization such 
that its parameterization metric has a positive entry in the right-upper corner. Therefore from 
CorollarylUit follows that (n-3, 1) $ NS(M(i)) . Since (n— 2, 2) -< (n— 2, 1) we can conclude from 
(|3~i2|) and (|3~i5)) that (n - 2, 2) is the minimal element of NS(M (t)) . In other words, (n - 2, 2) 
is the characteristic pair of our classification problem. Also, it implies that the characteristic 
matrix C of the problem is equal to Norm(P), where P is as in Theorem |^1 If n = 5 then 
obviously C = Norm(P) = P. For n > 5 one can calculate all nontrivial entries of Norm(P), 
using identities (|2. . It gives all nontrivial entries of the characteristic matrix C: 

C n —l,n-4: = n — 2, 
C n -2,n-A = n(n - 3), 

n— 3— i 

c 4 ,„_ 4 = n(i + £ rr) + (££?) + ; /) + (*-£?)). 2 < z < n - 3 , 



(3.44) 



!=1 

Cl, n -4 = n, 

cii = -(( n 7-t J ) + 2(750 + 2(7-*fl + (7!r) - 1- 

-Em)- l<J<n-5. 
1=1 

Recall that Cij is the intrinsic number of the functional invariants of i + 1 variables and the 
weight j + 2. In order to obtain a characteristic parameterization from the parameterization 
by the tuple of the primary invariants one can implement some series of rearrangement of the 
primary invariants according to the series of elementary transformations from the matrix P to 
Norm(P), as was described in the previous section (see, for example, the formula (|2.11j) and the 
paragraph after it). 

Remark 6 Two control systems y = J-(y,v) and y = J-(y,v), with m-dimensional state- 
space S and one dimensional control space V are called micro-locally state-feedback equivalent 
at the point (yo,vo) E S x V, if there exist the state-feedback transformation 

= $l(2/o) 5 V = $2(2/0, «o) 

such that in the neighborhood of (2/0,^0) in S x V the following identity holds: 

^(^GO^aG/.u)) =d$i.F(i/,tO. 

The afline (1, m + 1) control system (|1.5|) will be called i/ie affine extension of the control system 
p.4j) . It is not difficult to show that two control systems with scalar input are micro-locally 
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state- feedback equivalent at some point (yo, vo) € S x V if and only if their affine extensions 
are locally equivalent w.r.t. the state-feedback transformations of the type (|1.3[) at the same 
point. Note also that the affine extensions of generic m-dimensional control systems with scalar 
input are generic in the set of all (l,m+ l)-affine systems. Using this fact and Theorems ^ 
121 one obtains the micro- local parameterization of non-affine m-dimension control systems with 
scalar inputs and m > 3 by the tuples of the primary invariants of their affine extensions (in 
C°° category for m = 3 and C u category for m > 4). Obviously, the Poincare series, the 
characteristic pair and the characteristic matrix of the micro-local state-feedback classification 
problem for m-dimensional control systems with scalar input are exactly the same as in the 
case of the state-feedback classification problem of (l,m + 1) control-affine systems. Besides, 
since in our method of normalization of (1, n) control-affine system with n > 5 we rectify the 
vector field /i, then in the case m > 4 a generic m-dimensional real analytic control system 
with the prescribed tuple (|3.38|) of the primary invariants of its affine extension is micro-locally 
state-feedback equivalent in the real analytic category to the system 

y s = fs(yi,---,y m ,v), i<s<m, 

such that f s (x 2 ,x 3 , . ..x m ,xi) = a s+1 (xi, . ..x m+1 ), wh ere th e tuple {a k (xi, . . . ,x m )}^^ is the 
solution of the system of partial differential equations (|3.34|) with boundary conditions, which 
can be expressed by the primary invariants ip^, and 4>kjh using (|3.37l) . 

4 Reduction of control-affine systems with two-dimensional in- 
put to the scalar input case in dimensions four and five 

For r = 2 the system has the following form 

? = /o(«)+ttl/l(?)+U!/!(?) ) ?6M, Ul , U2&R. (4.1) 

Our aim is to assign to the system (|4.1j) in a canonical way an affine subsystem with scalar 
input 1 . It turns out that in the case n = 4 the original system can be recovered from it uniquely 
up to a feedback transformation, while in the case n = 5 such unique recovering is possible 
after introducing an additional invariant function of n variables (which is natural in view of the 
estimates for the number of functional parameters, given in the Introduction). 

4.1 Preliminaries. Let us look on (|4.1|) . as on the time optimal control problem and find 
its extremals. First we introduce some notations. Let T*M be the cotangent bundle of M with 
canonical symplectic form a. Denote by hi, < i < 2, the following functions on T*M: 

h{\) = P -fi(q),X={p,q),qeM,pe T*M. (4.2) 

For a given function G : T*M i— > K denote by G the corresponding Hamiltonian vector field 
defined by the relation cr(G, •) = —dG(-). For a given vector distribution D on M (i.e., a 
subbundle of the tangent bundle), define the Ith power D l by the recursive relation 

D l = D 1 - 1 + [D, D l -\ D 1 = D, 

and denote by (D ) <ZT*M the annihilator of D\ namely 

(D 1 ) 1 = {(p, q) G T*M :p-u = 0VuG D l (q)}. 

the meaning of the word "subsystem' is that at any point q the set of its admissible velocities is a subset of 
the set of the admissible velocities of the original system at q. 
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In the introduced notations the Hamiltonian of Pontryagin Maximum Principle for the time 
optimal problem (|4.1|) can be written as follows: 

H(X,u 1 ,u 2 ) = h (X)+u 1 h 1 (X)+u 2 h 2 (X), A G T*M, Ul ,u 2 G R. (4.3) 

Let 7(-) be an extremal of (|4,1[) with extremal control functions ui(t) and u 2 (t). Then 

7(t) = h {j(t)) +u 1 (t)hi{"/(t)) +Mt)H"f( t )) ( 4 - 4 ) 
and from the maximality condition for H it follows that 

7 (-) C {A G T*M : h^X) = h 2 (X) = 0}. (4.5) 

If we denote D 2 = span(/i,/2), then (|4.5j) is equivalent to 7(-) C (Z^) -1 - Combining (|4.4j) and 
(|4.5|) . we obtain 

& r(f) h i [i(t))=0,i = l,2. (4.6) 
Then from (|4.4j) and (|4.6() it follows 



{/to, MW*)) + «i(t){fci,/i2}(7(t)) = 



(4.7) 



(here are Poisson brackets of the Hamiltonians /ij and hf. {h{,hj} = dhj{hi)). Now 

suppose that 

dim D\ = 3 (4.8) 

Then relations (|4.7j) implies that the extremals of (|4.1|) . lying in (Z^^V-Dl)" 1 ", are exactly the 
integral curves of the vector field 

{/i 2 ,/ii} {/i2,ni} 

{/to, /12} 

(which is the Hamiltonian vector field, corresponding to the Hamiltonian X = /jq + 77 — ; — 7/11 + 

{h 2 ,hij 

{h±, ho} rm 

— — 7—rh 2 ). Denote by V the affine subbundle of TM, defined by system (|4.1j) and V(q) be the 

{^2,^1} 

set of all admissible velocities of the system (|4.1jl at the point q, 

V(q) = {fo(q) + + u 2 f 2 {q) : u u u 2 G M}. 

Let 7r : T*M wMbe the canonical projection. The set 

Ext(g) = {tt*X(\) : A G T*M D (Z^VZA^}, g G M (4.10) 

is the subset of V(q), consisting of the velocities of all extremal trajectories starting at q and 
having a lift in (D^^D^ . 

Among all extremals on (D 2 ) ± \(D2)~ L , one can distinguish so-called abnormal extremals, i.e., 
the extremals lying on the zero level set of the Hamiltonian X. Denote D% = span(/o, fi, f 2 ) 
and suppose that 

dim (Dl + D a ) = 4 (4.11) 

The set 

Abn(g) = {vr*X(A) : A G T*M n ( J D 3 ) ± \Pi) ± }, q € M (4.12) 
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is the subset of Ext(g), consisting of the velocities of all abnormal extremal trajectories starting 
at q and having a lift in (D 2 )- L \(D2) ± ■ One can show that for generic affine systems of the type 
(|4,1|) Ext(g) = V(q) in the case n > 5 and Abn(g) = V(q) in the case n > 6. But in the case 
n = 4 and n = 5 either Ext(g) or Abia(q) (or both of them) define the proper subsystem of the 
original system (|4.1|) . Moreover, it turns out that these subsystems are affine with scalar input, 
so one can apply the theory of the previous section. Now let us consider the cases n = 4 and 
n = 5 separately. 

4.2 The case n = 4. Let 



[V,D 2 ](q) = {[X,Y](q) : X £ V, Y G D 2 , are vector fields}. 



It is not difficult to show that [V, D 2 ](q) is a linear space and 

[V J D 2 ](q)= S p a n(f 1 (q)J 2 (q),[f u f 2 ](q),[f J 1 ](q),[f ,f 2 ](q)) (4.13) 
The crucial observation is formulated in the following 



Proposition 3 The set Ext(g) is an affine line, provided that j4.8j) holds and 

dim[V, D 2 ](q) =4 (4.14) 



Proof. Take some vector field such that the tuple (/o, /i, /b, /j) constitutes the frame on 
M. Denote by the structural functions of this frame, i.e., the functions, satisfying 

3 

[/i,/;] = E4/fc- ( 4 - 15 ) 

k=0 

Using the following well-known property of the Poison brackets 



{hi, hj}(p, q)=p- [fi, fj}(q), qeM,pe T*M 
and (|4.9|) . one can easily obtain that 



(4.16) 



Ext(q)=U(q)nV(q), (4.17) 

where 

n(g) = {{4 2 {q)v + c 3 12 (q)v)Uq) + (c° (^+ 

4)(?)/')/i(?) + ('(!i(?) I ' + ^i(?)/i)/j(?):/i ) i'eR} 

From assumption ()4.14|l and identity (|4.13|) it follows that H(q) is a plane. Assumption (|4.8|) 
implies that the plane n(g) is not parallel to the plane V(g). Note also both H(q) and V(g) 
belong to Ds(q). Hence by (|4.17|) the set Ext(g) is an affine line. □ 



Consider the control system such that Ext (q) is its set of the admissible velocities at q. 
By Proposition it is an affine system with scalar input. We call this system the reduction 
of the four- dimensional control- affine system (14.11) . The following proposition gives another 
characterization of the reduction of the system (|4.1(l : 
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Proposition 4 Assume that the four- dimensional control- affine system \4-l\) satisfies the 
conditions \4-$[l and \4-14\l - Then the subsystem 

q = go + ug 1 , (4.19) 

of \4-l\ ) * s its reduction if and only if 

[go,gi]£D 2 . (4.20) 

Proof. By definition, the system 1)4.19)1 is the reduction of 1)4.1)1 if and only if 

Ext(g) = {g {q) + t gi (q) : t E R}. (4.21) 

On the other hand, on can take from the beginning /p = go and fx = g\. Then comparing 1)4.21)1 
with (|4.17f) and ()4.18)1 we obtain that the system (|4.19f) is the reduction of l)4.1|l if and only if 
c oi = c oi = 0' which is equivalent to [<?p,<?i] £ s P an (<?i, /2) = D 2 . □ 

Corollary 2 Assume that a four- dimensional affine control system j4.iyQ satisfies 

dim span (51, [gi,g ], [gi, [gx, g ]], [g , [gi,go]]) = 4 (4.22) 

Then any four- dimensional control- affine system with two-dimensional input, having the system 
14-iyQ as its reduction, is feedback equivalent to the system 

q = go + uxgi + u 2 [g ,gx], Ux,u 2 eR. (4.23) 

Proof. First by assumption 1)4.22)1 and relation 1)4.13)1 the system 1)4.23)) satisfies conditions 
()4.8|) and 1)4.14)1 (where /o, fx, and f 2 are replaced by go, g\ and [<?o> <7i])- Hence, by Proposition 
01 the system 1)4.23)1 admits the reduction and by Proposition 21 this reduction is the system 
()4.19)1 . On the other hand, suppose that some system 1)4.1)1 has the reduction ()4.19)) . Then from 
the previous proposition [<?o>Si] £ span(si) f 2 ). According to ()4.22)) . g\ and [5o>ffi] are linearly 
independent. Hence the system (|4.1j) is feedback equivalent to 1)4.23)) . □ 

According to the previous proposition, a generic four-dimensional control-affine system with 
two-dimensional input can be uniquely, up to a feedback transformation, recovered from its 
reduction. Suppose that the reduction (|4.19|) of the system ()4.1I) satisfies (|4.22f) and 

dim span (g , gi, [gi, go], [gi, [01,50]] )= 4. (4.24) 

Then we can apply to the system 1)4.19)1 all constructions of section 2. In particular, one can 
construct the tuple of the primary invariants of 1)4.19)1 at a given point, which are also feedback 
invariants of the original system (|4.1j) . Note that the set of germs of systems of the type 
1)4.1)1 having the reductions, which satisfies conditions 1)4.22)1 and 1)4.24)1 . is generic. Combining 
Theorem^ Corollary [2 and normal form ()3.32|) . we obtain the following classification of generic 
germs of systems of the type 1)4.1)1 in terms of the tuple of the primary invariants of their 
reductions: 

Theorem 3 Given two arbitrary germs I\ and X 2 of functions of four variables at 0, four 
arbitrary germs (3\, f3 2 , /?3, /?4 of functions of three variables at 0, and four arbitrary germs 
■0i> ^2 j ^4 of functions of two variables at there exists a unique, up to state-feedback transfor- 
mation of the type 11. fy) . four- dimensional control-affine system with two-dimensional input such 
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that its reduction satisfies genericity assumptions &4.2Hfy and M-24\ ) an d (Ti,l2, 0i, f3 2 , @3, Pi, i>i, 
^2, ip3, ^4) is the tuple of the primary invariants of the reduction at the given point go- This 
control system is state-feedback equivalent to the following one: 

x 1 = l + (N + B 1 p 2 )u 1 + + A§jr)«2 

x 2 = (Pi + B 2P2 ) Ul + ^ +J B 2 |^)u 2 , iii,« 2 6l, (4.25) 
Xi = Bip 2 ui + [3i^P L u 2 , i = 3,4 

OX\ 

where N is the solution of hH.2!J\) . pi, i = 1,2 are the solutions of i'-i.S(J\) . and B^, 1 < k < 4, 
are as in The Poincare series, the characteristic pair and the characteristic matrix of 

the classification problem are exactly the same as in the case of (1,4) control- affine systems. 

Remark 7 It is easy to show that in the case n = 4 the set Abn(g) consists of one vector 
provided that (|4.11l) holds. Besides, if the system (|4.19f) is the reduction of the system (|4.1j) and 
it satisfies (|4.22jl and (|4.24|) . then Abn(g) is exactly its canonical drift. □ 

Remark 8 Actually, there is another intrinsic way to assign to the system (|4.1|) . satisfying 
(|4.11|) . an affine subsystem with scalar input: As a drift one can take again Ab(q). It remains 
to define canonically the direction of the affine line of the reduction. For this note first that 
the distribution D 2 satisfies (|4.8j) because of assumption (|4.11|) . Therefore through any point 
of M the unique (unparameterized) abnormal extremal trajectory of the rank 2 distribution 
D 2 passes: the line subdistribution L of D 2 , tangent to the abnormal extremal trajectories at 
any point is characterized by the relation [L,D 2 ] C D 2 . The direction of the affine line of the 
reduction can be taken parallel to L. The direction of L is different in general from the direction 
of the affine line in the first reduction. But this new reduction is worse than the previous one, 
because the original system (|4.1[) is not uniquely recovered from it: if (g~o,g~i) is the canonical 
pair of the new reduction (by construction and the previous remark go = Abn), then the field 
f 2 can be taken in the form f 2 = ago + [<h, <7o]> where the function a satisfies some second order 
ordinary differential equation along each integral curve of g\. Note that the direction g\ depends 
on the second jet of the original system ()4.1|) . while the direction of the affine line in the first 
reduction depends only on the first jet. This could be the reason for the loss of some information 
about the original system during the reduction described in the present remark. □ 

4.3 The case n = 5. In this case by analogy with Proposition |31 we have 

Proposition 5 The set Abn(g) is an affine line, provided that J^. i j| ) holds and 

dim Dl = 5. (4.26) 

Proof. Take some vector fields ^3 and f^ such that the tuple (/o ; fi, f2, f3, fi) constitutes 
the frame on M. By analogy with (|4.15|) . let c^, < i,j,k < 4, be the structural functions 
of this frame. From Q4.9j) and (|4.12|) . using (|4.16|) and the fact that ho = on (D^)^, one can 
easily obtain that 

Abn(q)=n 1 (q)nV(q), (4.27) 

where 

Ili(g) = {{c\ 2 {q)v + c\ 2 {q)p)h{q) + {c 3 20 (q)v+ 

4 (q)p)fi(q) + (4i(q> + 4l(<l)v)f2(q) :ftf£R) 
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From assumption Q4.26JI and identity Q4.13JI it follows that IIi (q) is a plane. Assumption (|4.11j) 
implies that the plane H-i(q) is not parallel to the plane V(q). Note also both Hi(q) and V(q) 
belong to D 3 (q). Hence by 1)4.27)1 the set Abn(g) is an affine line. □ 

Consider the control system such that Abn(g) is its set of the admissible velocities at q. 
By Proposition |S] it is an affine system with scalar input. We call this system the reduction 
of the five- dimensional control- affine system 1)4. 1)1 . The following proposition gives another 
characterization of the reduction of the system (|4.1|) : 

Proposition 6 Assume that the five- dimensional control- affine system \4-.l\j satisfies the 
conditions \4-2b\j and \4-H\j - Then the subsystem 

q = go + ugi, (4.29) 

of \4-l\) * s its reduction if and only if 

\ga,gi]eB 3 . (4.30) 

Proof. By definition, the system 1)4.19)1 is the reduction of 1)4.1)1 if and only if 

Abn(qf) = {g {q) + t gi {q) : t G K}. (4.31) 

On the other hand, on can take from the beginning fo = go and f\ = g\. Then comparing 1)4. 31|) 
with (|4.27j) and ()4.28)) we obtain that the system 1)4.29)1 is the reduction of 1)4.1)1 if and only if 
c oi = c oi = °, which is equivalent to [50,51] G span(# , 5i, ^2) = D 3 . □ 

Corollary 3 Assume that a five- dimensional affine control system i4.19\ ) satisfies 

dim span (50, 51, [5i,5o], [51, \gi,9o\], [go, [51, 5o]]) =5 (4.32) 

Then a five- dimensional control- affine system with two-dimensional input, has the system \4-iyj) 
as its reduction if and only if it is feedback equivalent to the system 

q = go + u 1 g 1 +u 2 (ag + [go,gi}), ui,u 2 £R, (4.33) 

where a is some function. 

Proof. First by assumption 1)4.32)1 the system 1)4.33)1 satisfies conditions 1)4.26)1 and 1)4.1 1)1 . 
Hence, by Proposition |S] the system ()4.33)1 admits the reduction and by Proposition this 
reduction is the system 1)4.19)1 . On the other hand, suppose that some system 1)4.1)1 has the 
reduction (|4.19|) . Then from the previous proposition [50,51] G span(c/o, 51, f 2 )- According to 
1)4.22)1 . 50, 51 and [50,51] are linearly independent. Hence, f 2 = £o5o + 651 + 6[5o,5i], where 
^2 / 0. Hence by a feedback transformation we can replace f 2 with ago + [50,51]- D 

So, in contrast to the case n = 4, five-dimensional control-affine systems with two-dimensional 
input cannot be recovered from its reduction only. Suppose that the system 1)4.1)1 has the re- 
duction 1)4.29)1 satisfying condition 1)4.32)1 and also 

dim span (50, 5i ! [5i,5o], [51, [51, go]), [51, [51, [5i,5o]]])= 5, (4.34) 
dim span(gi, [g , [50, 5i]] , [51, 5o], [51, [5i,5o]], [51, [fll, [5i,5o]]]) = 5. (4.35) 
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Then we can apply to the system (|4.29[) all constructions of section 2. In particular, let (Go, Gi) 
be the canonical pair of the system (|4.29|) . As before, let V be the affine subbundle of TM, 
defined by system (|4.1|h Then by the same arguments, as in the proof of Corollary 01 there exist 
a unique vector field G2 € V such that 



G2 — RGq + [Gq, G{\. 



(4.36) 



By construction, the function R is a feedback invariant of system ()4.1|) . Moreover, the system 
(|4,ll) can be uniquely, up to a feedback transformation, recovered from its reduction and the 
function R. 

Fix some point qo in M. Let $5 be as in (|3.17|) . Denote 



TZ = R o $ F 



(4.37) 



TZ is the germ of a function of five variables at 0. We call it the recovering invariant of (|4.1|) at the 
point q$. By construction, it is invariant of the weight 0. Note that the set of germs of systems 
of the type (|4,1|) having the reductions, which satisfies conditions ()4.32j) . 1)4,34}) . and ()4.35j) . is 
generic. Using Theorem |2] in the case n = 5 and the definition of the recovering invariant, we 
obtain the following classification of generic real analytic germs of systems of the type (|4.1j) in 
terms of the tuple of the primary invariants of their reductions and their recovering invariant: 



Theorem 4 Given 4 germs T\, I2, I3 



and TZ of real analytic functions of 5 variables at 
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such that 1^(0) 7^ 0, 10 germs {/3fc}| =1 an d { < Pk22}k=i °f rea ^ analytic functions of A variables at 
0, 10 germs {ipk}k=i and {4>k2z]\ = i of real analytic functions of 3 variables at 0, and 5 germs 
{<fik24}k=i of real analytic functions of 2 variables at there exists a unique, up to state-feedback 
real analytic transformation of the type kl.ty . five- dimensional real analytic control- affine system 
with two input such that first its reduction satisfies genericity assumptions [4.. Sty , \4-^4\) ; \4-35\) ; 
secondly ({Zj}j = i, {Pk}\ = \-i {^k}\ = i, {<Pk2l : 1 < k < 5,2 < I < 4}^ is its tuple of the primary 
invariants, and finally TZ is its recovering invariant at the given point q$. In other words, 
the tuples of the primary invariants give the regular (3, 2) -parameterization of the considered 
classification problem in the real analytic category with the following 4x4 parameterization 
matrix P: 



P 



The Poincare series M(t) of the considered classification problem satisfies 



(° 








5 \ 











10 











10 


\1 








3/ 



(4.38) 



M(t) 



1 



(i-ty 



■+ 



10 



10 



■+ 



(l_ t )S (l_ t )4 (1_£)5 (1 _ t) 



(4.39) 



By the same arguments, as in the case of (l,re)-affine control system with n > 5, treated 
in section 3, one can show that (3,2) is the characteristic pair of the considered classification 
problem (just use formulas (|2.1fij) . (|2.17j) . and Corollary [J). The characteristic matrix C of the 
problem is equal to Norm(P), which can be found easily by the series of elementary transfor- 
mations. Namely, 



C = Norm(P) 



A 


3 


6 


5 \ 











16 











13 


\o 











(4.40) 
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Conclusion 2 The characteristic parameterization of (2, 5) control- affine systems, up to 
state-feedback transformations, consists of 4 functional invariants of 5 variables and the weight 
3, 13 functional invariants of A variables and the weight 3, 16 functional invariants of 3 variables 
and the weight 3, 1 functional invariants of 2 variables and the weight 0, 3 functional invariants 
of 2 variables and the weight 1, 6 functional invariants of 2 variables and the weight 2, and 5 
functional invariants of 2 variables and the weight 3. 

In order to obtain a characteristic parameterization from the parameterization by the tuple 
of the primary invariants and the recovering invariant one can implement some series of rear- 
rangement of the primary invariants according to the series of elementary transformations from 
the matrix P to Norm(P), as was described in the previous section (see, for example, formula 
(|2.11j) and the paragraph after it). 
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